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FOREWORD 


This report has been prepared to expedite early domestic dissemination of 
the information generated under the contract. The data and conclusions must 
be considered preliminary and subject to change as further progress is made on 
this program. This is a progress report covering the work done during the 
second 12 months of the contract, and it is not a final report. 




TABLE OF CONTENTS 


Section Page 

1.0 INTRODUCTION 1 

2.0 TECHNICAL PROGRESS 7 

2.1 Task I - Literature Survey 7 

2.2 Task II - Design of Structural Analysis 

Software Architecture 12 

2.3 Task III - Thermodynamic Engine Model 12 

2.4 Task IV - Software Development 

2.5 Task V - Mission Model Development 12 

2.5.1 Component Temperature and Pressure Decompo- 
sition for Hot Synthesis 12 

2.5.2 Turbine Blade and Vane Temperature and 

Pressure 17 

2.6 Task II - Structural Analysis Methods Evaluation 51 

2.7 Task VII - Thermodynamic Loads Model 52 

2.8 Task VIII - Component Specific Model Development 60 

2.8.1 Geometric Modeling 60 

2.8.2 Remeshing and Mesh Refinement 66 * 

2.8.3 Self-Adaptive Solution Strategies 82 

APPENDIX A - TASK II - DESIGN OF STRUCTURAL ANALYSIS 

SOFTWARE ARCHITECTURE 96 

APPENDIX B - COMPONENT TEMPERATURE AND PRESSURE 

DECOMPOSITION AND SYNTHESIS PLAN 107 


v 



LIST OF ILLUSTRATIONS 


Figure Page 

1. Component Specific Modeling Base Program. 4 

2. Component Specific Thermomechanical Load Mission Modeling. 4 

3. Component Specific Structural Modeling. 5 

4. Material Thickness Temperature Gradient. 15 

5. Coordinate System for Cooling Effectiveness. 16 

6. HPT Stage 1 Blade Cooling Effectiveness. 18 

7. F102 Bl-B HPT Blade Cooling Effectiveness. 19 

8. Load Cycles 1, 2, 3, and 4. 26 

9. Shingle Segment. 30 

10. CYANIDE Model. 30 

11. Temperature Distribution on Shingle at Peak Condition. 31 

12. Nonlinear Patch in Elastic Matrix. 37 

13. Data Points Used for Interpolation Routines in the TEM. 55 

14. Quadrant Error Distribution. 56 

15. Typical Nugget. 61 

16. Combustor Recipe Process. 62 

17. Combustor Nugget Element Models. 63 

18. Cooled Turbine Blade Cross Section. 64 

19. Plot of Exerior and Interior Points for a Typical Airfoil. 65 

.20. Points Chosen for a Spar Configuration. 67 

21. Coarse Mesh Selection for a Blade With Cavities. 68 

22. Coarse Mesh Selection for a Blade Without Cavities. ' 69 

23. Recipe Parameters and Master Region Model for Plate With 

Hole. 74 

24. Generated Master Region. 75 

25. Coarse Discretized Model. 76 

26. Magnification of Hole Region for Coarse Model. 77 

27. Refined Discretized Model. 78 

28. Magnification of Hole. Region of Refined Model. 79 

vi 



LIST OF ILLUSTRATIONS (Concluded) 


Figure Page 

29. Coarse Mesh. 80 

30. Medium Mesh. 80 

31. Fine Mesh. 80 

32. Effect of Mesh Density on Stress Prediction. 81 

33. Problem Boundary Conditions and FEM Meshes. 83 

34. Results of the Effect of Mesh Density. 84 

35. Cycle Crack Growth Model. 89 


vii 



LIST OF TABLES 


Table Page 


I. Linear Fit Constants for Equation 8. 14 

II. Data Base. 24 

III. Plastic Strain as a Function of Temperature 

and Pressure. 25 

IV. Details for the Four Load Cycles. 27 

V. Predicted Plastic Strain Response for the Fan Load Cycles. 27 

VI. Comparison of the Plastic Strain Prediction Methods. 28 

VII. Case 1 Results - R Components of Strain. 33 

VIII. Case 3 Results - R Components of Strain. 34 

IX. Case 2 Results - R Components of Strain. 35 

X. Case 1 Results - Z Components of Strain. 36 

XI. Case 2 Results - Z Components of Strain. 37 

XII. Case 3 Results - Z Components of Strain. 38 

XIII. Turbine Blade Tip Model. 40 

XIV. Case 1 Results - Z Components of Strain. 41 

XV. Case 1 Results. 44 

XVI. Case 2 Results. 45 

XVII. Total Strain. 46 

XVIII. Creep Strain. 46 

XIX. Eight-Element Patch Total Strain. 47 

XX. Eight-Element Patch Creep Strain. 47 

XXI. Error Summary. 57 

XXII. CF6-50C2 Validation Case Error Analysis. 58 

XXIII. Stability of Dynamic Time Incrementing. 87 

XXIV. Displacement Tolerance Criteria. 89 

XXV. Global Euclidian Criteria. 90 

XXVI. Percentage Change Criteria. 90 

XXVII. Utilization of L(^) Vector Norms. 94 

XXVIII. Utilization of 1 ( 2 ) Vector Norms. 94 

XXIX. Utilization of LC^) Vector Norms. 95 


viii 



1.0 INTRODUCTION 


Modern jet engine design imposes extremely high loadings and temperatures 
on hot section components. Fuel costs dictate that minimum weight components 
be used wherever possible. In order to satisfy these two criteria, designers 
are turning toward improved materials and innovative designs. Along with 
these approaches, they also must have more accurate, more economical, and more 
comprehensive analytical methods. 

Numerous analytical methods are available that can, in principle, handle 
any problem that might arise. However, the time and expense required to pro- 
duce acceptable solutions is often excessive. This program addresses this 
problem by setting out a plan to create specialized software packages which 
will provide the necessary answers in an efficient, user-oriented, stream- 
lined fashion. Separate component-specific models will be created for burner 
liners, turbine blades, and turbine vanes using fundamental data from many 
technical areas. The methods developed will be simple to execute, but they 
will not be simple in concept. The problem is extremely complex and only by a 
thorough understanding of the details can the important technical approaches 
be extracted. The packaging of these interdisciplinary approaches into a 
total system must conform to the modular requirements for useful computer pro- 
grams . 

The overall objective of this program is to develop and verify a series 
of interdisciplinary modeling and analysis techniques that have been special- 
ized to address three specific hot section components. These techniques will 
incorporate data as well as theoretical methods from many diverse areas 
including cycle and performance analysis, heat transfer analysis, linear and 
nonlinear stress analysis, and mission analysis. Building on the proven tech- 
niques already available in these fields, the new methods developed through 
this contract will be integrated to provide an accurate, efficient, and uni- 
fied approach to analyzing combustor burner liners, hollow air-cooled turbine 
blades, and air-cooled turbine vanes. For these components, the methods 
developed will predict temperature, deformation, stress, and strain histories 
throughout a complete flight mission. 
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This program, to a great extent, draws on prior experience. This base of 
experience is invaluable for understanding the highly complex intersctions 
among all the different technical disciplines as well as for estimating the 
importance of different engine parameters. In particular, there are four spe- 
cific areas in which experience is especially beneficial. 

First, with the recent increases in fuel costs, greater emphasis has been 
placed on* more accurate solutions for stresses and strains in order to under- 
stand and improve the durability and life of hot section components; Conven- 
tional linear elastic analyses are no longer sufficient instead, they now 
provide the boundary values for more refined creep and plasticity calculations. 
These nonlinear analyses are now performed routinely as part of the design 
process at General Electric. This extensive experience with these plasticity 
and creep methods contributes directly to developing component specific models. 

Second, advances in 3-D modeling capability are being achieved by the 
concepts developed under the NASA-supported ESMOSS program. ESMOSS concepts 
provide the basis to develop an efficient modeling system for geometric and 
discretized models of engine components. 

Third, the NASA- funded Burner Liner The rmal/ Structural Load Modeling Pro- 
gram contributes strong support to this program. The specific area addressed, 
transfer of data from a 3-D heat transfer analysis model to a 3-D stress analy- 
sis model, will provide the background and framework for the data interpolation 
required for all thermomechanical models in this contract. 

Fourth, over the past 10 years General Electric has developed internally 
a family of computer programs: LASTS, OPSEV, and HOTSAM; These programs all 

have the common thread of using selected points from cycle data, heat transfer, 
and stress analyses and a decomposition/ synthesis approach to produce accurate 
values of temperature, stress, and strain throughout a mission. These programs 
are totally consistent with the overall objectives of this program, and repre- 
sent a proven technology base upon which the component specific models are 
being developed. Significant advances being made are the inclusion of non- 
linear effects and the introduction of improved modeling and data transfer 
techniques . 
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The program is organized into nine tasks which can logically be separated 
into two broadly parallel activities (Figure 1). On the right of Figure 1 we 
have the Component Specific Thermomechanical Load Mission Modeling path. 

Along this path a Decomposition/Synthesis approach is being taken. In broad 
terms, methods are being developed to generate approximate numerical models 
for the engine cycle and the aerodynamic and heat transfer analyses needed to 
provide the input conditions for hot parts stress and life analysis. 

The left path, Component Specific Structural Modeling provides the tools 
to develop and analyze finite element nonlinear stress analysis models of com- 
bustor liners and turbine blades and vanes. These two paths are shown in more 
detail in Figures 2 and 3. 

Software Development, Task IV, consists of planning and writing the com- 
puter programs for both paths, with the necessary interconnections, using a 
structured, top down approach. 

In the Thermomechanical Load Mission Modeling portion of the program 
(Figure 2) , we are developing in Task III a Thermodynamic Engine Model which 
generates the engine internal flow variables for any point on the operating 
mission. The method for doing this is described below. Task V is developing 
techniques to decompose flight missions into characteristic mission segments . 
In Task VII a Thermomechanical Mission Model is being developed. This uses 
the flow variables from the Thermodynamic Model to determine metal temperature 
and pressure distributions for a representative combustor liner and turbine 
blade and vane. 

Individual tasks for the Structural Modeling activity are shown in Figure 
3. The requirements of Software Design, Task II, have been factored into Task 
VI, the evaluation of the structural analysis methods which were selected for 
evaluation in Task I. Task VIII provides the capability for structurally 
modeling current state-of-the-art combustor liners and hollow turbine blades 
and vanes, given the defining dimensional parameters. These parameters will 
be chosen to facilitate parametric studies. 

The component specific models are being developed in two steps . In the 
first a geometric modell is defined. In the application of the Component 
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Figure 1. Component Specific Modeling Base Program, 



Figure 2. Component Specific Thermomechanical Load 
Mission Modeling. 
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Figure 3. Component Specific Structural Modeling. 
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Specific Modeling Program these data are then transferred to the Thermomechan- 
ical Load Mission Model to provide the geometry for determining component 
pressures and temperatures. Thus, a data transfer link is being developed to 
do this in Task IV, Software Development. The capability for generating from 
the geometric model a discretized, finite element model is also a part of Task 
VIII. At this point another link between the two paths is needed to transfer 
the component temperatures and pressures from the Thermomechanical Load Model 
to the finite element model, interpolating the data as needed to define nodal 
temperatures and pressures. This also is being completed in Task IV. 
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2.0 TECHNICAL PROGRESS 


2.1 TASK I - LITERATURE SURVEY 

The first task of this program was to perform a literature survey of 
available methods, techniques, and solution strategies that can be used to 
geometrically model, display, and structurally analyze burner liner, turbine 
blades and vanes. NTIS, NASA, DTIC, and internal General Electric Company 
documents were searched. As a result of this survey, 85 papers and 8 books 
and procedures were discovered with pertinent information. As a result of 
evaluating this information, recommendations were made on the technology to be 
incorporated into the base and advanced portions of this program to the NASA 
Program Manager. The NASA Program Manager approved the following program 
content: 

GEOMETRIC SHAPE GEOMETRIES 
Base Program 

• ESMOSS - considered commercial and university codes 

Advanced Program 

• ESMOSS 

LINEAR AND NONLINEAR FINITE ELEMENT METHODS 
Base Program 

• 8-, 16- , and 20-noded isoparametric elements - considered beams, 
plates, and shells 

Advanced Program 

• Supplements from 3D inelastic program 

• Decomposition and synthesis methods 
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SPECIALTY ELEMENTS 


Base Program 

• None - considered slave, large FEM, hybrid, embedded hole 


Advanced Program 

• Supplements from 3D inelastic program 

• Decomposition and synthesis methods 


NON-FINITE ELEMENT METHODS 
Base Program 

• • None - considered finite difference, boundary integral 
Advanced Program 

• Supplements from 3D inelastic program 

• Decomposition and synthesis methods 

STIFFNESS AND MASS MATRIX ASSEMBLY TECHNIQUES 
Base Program 

• Element level algorithm 


Advanced Program 

• Frontal solution methods 

• Supplements from 3D inelastic program 

MATERIAL BEHAVIOR CHARACTERISTICS 


Base Program Elastic 

• Isotropic X 

• Orthotropic X 

• Temperature-dependent X 

• Cyclic Plasticity 

• Creep 


Inelastic 

X 

X 

X 

X 
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Advanced Program 


Elastic 


Inelastic 


• Viscoplasticity X 

• Anisotropic X X 

SOLUTION TECHNIQUE 
Base Program 

• SESOL (sparse matrix solver, skyline storage scheme) 


Advanced Program 


• Vectorized COLSOL (active column solver, compacted storage) 

• Frontal solution methods 

LINEAR AND NONLINEAR SUBSTRUCTURES 
Base Program 

• Superelements 

• Conventional finite element 

Advanced Program 

• Multilevel sub structuring 

SOLUTION STRATEGIES 
Base Program 

• Dynamic time incrementation 

• Dynamic load incrementation 

• Solution acceleration schemes - Aitkens extrapolation. Over- 
relaxation scheme 

Advanced Program 

• Supplements from 3D inelastic program 

• Padovan's techniques 



CENTRIFUGAL STIFFNESS MATRIX FORMULATION 


Base Program 

• Two-step method 

Advanced Program 

• Large deformation method (updated Lagrangian) 

EIGENVALUE AND EIGENVECTOR EXTRACTION 
Base Program 

• Subspace iteration 

• Master nodes 

Advanced Program 

• Lanczos method 

AUTOMATIC DISCRETIZATION PROCEDURE 
Base Program 

• ESMOSS 

. • Master regions - considered commercial and university codes 

Advanced Program 

• ESMOSS 

REMESHING AND GRID OPTIMIZATION 
Base Program 

• Progressive subdivision 

• Constraint equations 

Advanced Program 

• Total realignment 

• Element upgrade 
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DATA DISPLAY 


Base Program 

• ESMOSS 

Advanced Program 

• NASA in-house 

PROVEN DATA SETS 

Stress-Strain Decomposition and Synthesis Techniques 

• Turbine blade nonlinear structural life analysis 

• Multiaxial cyclic thermoplasticity analysis with Besseling's sub- 

volume method 

Component Temperature and Pressure Synthesis Techniques 

• CF6-50 rolled ring combustor 

• CF6-80A high pressure sector test 

• CF6-6 single shank turbine blade 

• CF6-80 and 80A turbine vanes 

• Thermodynamic engine model 

• Thermomechanical loads model - build on in-house expertise: OPSEV, 

OPSEV-A, LASTS, HOTSAM 

Survey of 3-D Heat Transfer Codes 

• Used survey results from "burner liner thermal structural load 
modeling" 

• Use internal THT-D program for development work 

• Thermal loads transfer module uses MARC or SINDA available to NASA 
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2.2 TASK II - DESIGN OF STRUCTURAL ANALYSIS SOFTWARE ARCHITECTURE 


The software architecture was designed using the methodology developed on 
the ESMOSS program. This development was carried out by a team whose members 
provided expertise in all of the pertinent areas. The architecture approved 
by the NASA Program Manager is contained in Appendix A. In addition to the 
program architecture, the preprocessor and postprocessor attributes are 
defined. 

2.3 TASK III - THERMODYNAMIC ENGINE MODEL 

The Thermodynamic Engine Model (TDE) was completed and approved by the 
NASA Program Manager in 1983* . 

2.4 TASK IV - SOFTWARE DEVELOPMENT 

This task consists of planning and writing the computer codes for both 
paths of this program with the necessary interconnectors. As such, it is a 
continuous and ongoing effort with the substance being covered under the other 
task headings . 

2.5 TASK V - MISSION MODEL DEVELOPMENT 

2.5.1 Component Temperature and Pressure Decomposition and Synthesis 

Based on our efforts in 1983 and the developments described below, a 
Component Temperature and Pressure Decomposition and Synthesis Plan was 
approved by the NASA Program Manager. This plan is outlined in Appendix B. 

2.5. 1.1 Combustor Liner Temperature and Pressure Decomposition 
and Synthesis 

An expression for the temperature gradient through the material thickness 
can be derived from the cooling effectiveness, the compressor discharge 
temperature, and the combustor exit temperature. The temperature gradient 
through the material can be calculated from 


*R.L. McKnight, "Component - Specific Modeling," Annual Status Report, NAS3- 
23687, 1983. 
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( 1 ) 


T . T = 2/a1 

a h c k 

where 

T„ = Hot side metal temperature, ° F 
n 

T^, = Cold side metal temperature, ° F 


Q/A 


Heat flux through material 


Btu 
hr- ft 2 


t = Material thickness, ft 


K = Metal conductivity 


Btu 

ft 2 -hr-° F 


The heat flux can be calculated from 

( 2 ) 


(3) 

where 

2 

h = Convection heat transfer coefficient Btu/hr-ft -° F 
c 

hr = Equivalent heat transfer coefficient for radiation to casing 
Btu/hr-ft 2 -® F 

^Liner = liner temperature, 0 F 

Substituting the heat flux expression into the gradient equation (1) 

gives 


Q/A = (h c + hr)(T c - T 3 ) 


or it is proportional to (T^ ^ 


Q/A = Ch c + hr)(T Unec - T 3 ) 


T - T 
A H C 


(h c + hr)t 


(T Iiner ' T 3 ) 


using the equation for cooling effectiveness, 


y„ = 


T - T 
4 Liner 


c T 4 -T 3 


(4) 


(5) 
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( 6 ) 


An equation for (T^ ner - T^) can be written as follows: 
(T liner ' V = (1 * T , + <*c ‘ »> T 3 

Substituting Equation 6 into the expression (4) gives 


T H~ T 3 


(1 - y c ) t 4 + (y c - 1) T 3 


a 


(h c + hr)t 


(7) 


The convection term, h c> varies with pressure and thus the gradient through 
the material thickness should be correlated with pressure. 

A THTD analysis was done at several pressure conditions and the calcu- 
lated temperature gradients were plotted vs P3 for several axial locations and 
the results are shown in Figure 4. The locations and coordinates are shown in 
Figure 5. As shown in the figure, the gradient data are correlated with 
pressure. The constants m and b in the equation 

T H " T C T H ' T C 

T tiner ‘ T 3 = « ' 7? T y ♦ (h c - 1) I3 = .^3 + b (8) 

are tabulated in Table I. 


Table I. Linear Fit Constants for Equation (8). 


Location 

X, inches 

m 

b 

i 

0.094 

12.3 x 10-5 

0.100 

2 

0.438 

14.1 x 10-5 

0.061 

3 

0.654 

9.0 x 10“5 

0.061 

4 

0.854 

10.7 x 10-5 

0.092 

5 

1.114 

28.1 x IQ -5 

0.168 
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Coordinate System for Cooling Effectiveness, 


t i 


Given the combustor exit temperature, T , the compressor discharge pres- 
sure, P^, and the compressor discharge temperature, T^, the temperature gra- 
dient through material thickness can be calculated from Equation 8 using the 
cooling effectiveness and the constants from Table I. 

2. 5. 1.2 Turbine Blade and Vane Temperature and Pressure 
Decomposition and Synthesis 

The spanwise distribution of overall local cooling effectiveness has been 
completed for the CF6-50 and B-1B NESIP (similar to CFM56-2) Stage 1 HP 
turbine blades; The results are compared in the attached Figure 6 which is 
based on THTD runs for 15, 50, and 77 % of blade span. The ratio of 

H c local^c 50°/ is unit y definition at the 50% span location. At the other 
two spans, the p ratio is identified for each of the sixteen points around 
the airfoil. The curves have been terminated at the locus or the average r| 
ratio for each. span. 

For the CF6-50 blade, this procedure defined a single curve for the pres- 
ure and suction surfaces. However, the B-1B blade is better represented by a 
two-branch curve at the 77% span (Figure 7) . We have reviewed the two 
blade designs for possible explanations of this characteristic. There is no 
obvious single cause. It is undoubtedly the combined result of configuration, 
coolant circuitry, the application of film cooling and variations in gas-side 
heat-transfer coefficients. 

It appears best to allow for incorporating separate curves for the pres- 
sure and suction surfaces, with freedom to input these curves for different 
blade designs. This is probably the thing to do for the B-1B blade, defining a 
separate curve for the pressure surface between Points 2 and 6. Points 7 and 
8 appear to be represented quite well by the curve for the suction surface. 
Using the suction surface curve for Points 2 through 6 could overpredict the 
temperatures by about 135° F at the 77% span. 

2.5.2 Stress-Strain Decomposition and Synthesis 

The decomposition and synthesis of stresses, strains, and deformations is 
technically the most challenging portion of this program. It requires 
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innovative methods to produce usable results for burner liners, turbine blades, 
and vanes. Thus, our goal under this task has been to compile a library of 
possible decomposition and synthesis techniques and to assess their validity. 
Among the techniques being considered are the following: 


• Assume that the structure remains totally elastic at all stress 
levels and do the decomposition and synthesis based on an elastic 
"pseudostress . " 

• Assume that the structure is deformation-controlled (strain range 
invariance). The first level of decomposition and synthesis would 
be based on deformations (total strains). A second level of 
synthesis could then introduce the effects of plasticity and creep 
by using the material response characteristics to partition the 
total strain into elastic, plastic, and creep components. 

• Assume that the structure is load-controlled (stress invariance). 
Decomposition and a first level of synthesis would be based on load 
terms reflecting the centrifugal loadings and the temperature and 
pressure distributions. A second level of synthesis could then 
introduce the effects of plasticity and creep by using the material 
response characteristics to determine the elastic, plastic, and 
creep strains that would be caused by the total load. 

• Use simplified nonlinear finite element modeling to decompose and 
synthesize the stresses, strains, and deformations in terms of the 
set of analyzed mission components. These simplified models could 
be either one 2D or 3D element or a nonlinear substructure. These 
models could use boundary conditons from the detailed analysis or 
they could be run as an intimate part of the detailed analysis. 

• Apply the method of superposition for the decomposition and 
sysnthesis of stresses, strains, and deformations. This method 
would be investigated based on the following hierarchy of calculated 
parameters : 

deformations 
strains 
- stresses 

We will determine to what degree these parameters can be decomposed 
and synthesized by superposing the results from the individual load- 
ing functions (temperature, pressure, rpm) . 

• Use linear and nonlinear interpolation of the results of a detailed 
analysis for decomposing and synthesizing stresses, strains, and 
deformations. The interpolating parameters would be second-level 
predicted temperatures, pressures, and rpm's. 

• Form look-up tables of deformations, stresses, and strains as 
functions of temperatures, pressures, and rpm ! s. These tables would 
then be used to decompose and synthesize the mission cycles. 


20 



• Finally, generate from test data an empirical model relating 

stresses, strains, and deformations to temperatures, pressures, and 
rpm's. With this model, mission cycles could be decomposed and new 
ones synthesized. 

We have developed the following general methodology: 

A. Know the component in question. 

B. Know the thermomechanical load cycle. 

C. Know the total response of the generic component through a detailed 
finite element analysis . 

1. The detailed analysis will encompass total thermomechanical 
loading conditions. 

2. From the detailed analysis: 

a. Can locate and categorize the critical location(s) in a 
specific component 

• User defined 

• Automated procedure 

b. Categorize the stress/strain response at the critical 
location 

• Uniaxial 

• Biaxial 

• Triaxial 

c. Identify the generic mode of load application 

• Tension 

• Compression 

• Thermal 

• Pressure 

• RPM (Inertia) 

• Strain (Displacement) Control 

• Load Control 

- Consider the effect of load split and the super- 
position of loads 

d. Material Behavior 

• Elastic 

• Elastic/Plastic 

• Elastic/Creep 

• Elastic/Plastic/Creep 

• Flow Rule and Yield Criteria Used 

• Creep Formulation Used 

- Rank and/or prioritize the ability to synthesize 
desired material response 
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Decomposition Techniques: 

A. Assess the methodology used to-date in our attempt to decompose the 
thermomechanical response history of baseline models. 

B. Serving as baseline data sets are two models previously validated 
under other NASA programs. 

1. "Turbine Blade Tip Model" NASA CR-165268 

• Uniaxial strain cycling conditions 

• Simulation of complex strain/temperature cycling 

conditions 

2. "Inelastic Shingled Combustor Hot Spot" NASA CR-2278 

• Biaxial strain cycling conditions 

• Simulation of a very complex thermal cycle 

C. Have used CYANIDE, our in-house nonlinear finite element code to 
analyze such models under select conditions 

1. Have economically and efficiently created a preliminary data 
base. 

2. Have resorted to use of the code because we have little or no 
relevant experimental data having direct application to our 
specific needs. 

3. Have yet to determine if the content of the data base encom- 
passes all the information we may require in the synthesis of 
stress, strains and deformations. 

D. Are presently evaluating the concept of degenerating thermomechani- 
cal response into component parts. 

1. For our situation, is the superposition of thermomechanical 
response a valid technique and does it suit our needs? 

• Can the total response be characterized as a sum, a 
function of the sum, or some function of the individual 
response? 

• To date we have concerned ourselves only with elemental 
results - strains in particular. Can the concept be 
extended to include displacements and/or stresses? 

Synthesis Techniques: 

A. For each specific component, construct an algorithm to compute fixed 
point mission-time profiles of "local" stresses, strains, and 
displacements. Essential to this effort is the assumption that we 
need only synthesize the response on the local level. 
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B. Will base the synthesis techniques on a defined set mission-segment 

component-station characteristics . 

• This may include the TML matrix and the characterization 
. procedure noted in general methodology. 

C. Once the response at the critical location has been characterized, 

the synthesis of local response can proceed along assorted paths. 

• Generate empirical relations between loads and deformation 
states. 

• Generate closed form solutions. 

• Perform a parallel path solution utilizing another analysis 
code. 

D. Assemble the parameters defining local response and make them 

available to the user. 

This year's effort in the stress/strain decomposition and synthesis tech- 
niques has focused on the development of methods to predict structural 
behavior in response to a defined set of loads. These methods are based on 
the premise that the structural response to an arbitrary set of loads can be 
accurately predicted, provided that the current loads are perceived as being 
small pertubations of the loads used in a previous analysis. We are examining 
the use of linear and nonlinear interpolation routines to use the information 
obtained from a prior analysis to synthesize current nonlinear strain response. 

To circumvent the complex behavior exhibited in a two or three-dimensional 
analysis, it was decided to first limit ourselves to simple uniaxial condi- 
tions. For our work we would examine a uniaxial specimen under the following 
conditions : 

1. No time-dependent material behavior 

• Total strain can be separated into elastic and plastic 
components 

2. Consider only small inelastic and plastic components 

• Range from 0.4% to 0.6% 

3. Consider a limited temperature range 

• Range from 600° F to 1200° F 

4. Will examine only simple load cycles 
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Our in-house nonlinear finite element code was used to analyze a uniaxial 
model under various loading conditions; The results from the analyses that 
constitute our data base are tabulated and presented in Table II. Equations 
that express the plastic strain as a function of load and temperature were 
defined and are presented in Table III. These equations and linear interpola- 
tion were used to predict plastic strain response for various load cycles. 


Table II. Data Base. 


T/P 

600 

700 

800 

900 

1000 

1100 

1200 

Total Str 

a ins 







110 

4381 

4480 

4580 

4683 

4787 

4892 

5000 

120 

4995 

5101 

5207 

5312 

5417 

5521 

5625 

130 

5630 

5735 

5839 

5943 

6057 

6187 

6320 

140 

6314 

6442 

6572 

6705 

6840 

6979 

7158 

Plastic S 

trains 







110 

430 

468 

505 

542 

579 

615 

650 

120 

685 

724 

761 

795 

827 

855 

880 

130 

961 

993 

1023 

1049 

1084 

1132 

1180 

140 

1296 

1336 

1385 

1435 

1485 

1534 

1622 


Having the data base in place and the relationship for plastic strain 
defined, assorted load cycles could be defined to test the quality (accuracy) 
of the strain prediction techniques. Four load cycles are presented in Figure 
8. Table IV illustrates the results obtained when the plastic strains are 
defined as being a function of stress. In Table V, predicted strains are 
those determined via our synthesis techniques and the computed strains are 
those strains obtained directly from our finite element code. Table VI illu- 
strates the results obtained for cycle No. 4 when four separate methods of 
strain prediction are compared. 

In our investigation of stress-strain decomposition and synthesis tech- 
niques we continued to use the combustor liner shingle segment model presented 
in Footnote 1. Noting that the thermal stress-strain response of the model is 
biaxial, our efforts are directed toward examining the nonlinear stress-strain 
behavior (or the history) of a defined critical location. 
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Table III. Plastic Strain as a Function of Temperature and Pressure 


Plastic Strain = f (Temperature) 


£p = AT D 

Pressure (ksi) 


A 

B 

110 


9.33347E-6 

0.5977902 

120 


6.66123E-5 

0.3643609 

130 


1.515E-4 

0.2868058 

140 


1. 633E-4 

0.3208346 

Plastic 

Strain = (Pressure) 




£p = AP B 


Temperature (° 

11 

A 

B 

600 


2.59397E13 

4.522953 

700 


7.35291E-13 

4.317914 

800 


1.18710E-12 

4.140312 

900 


4.07319E-12 

3.982063 

1000 


7.94597E-12 

3.85249 

1100 


1.29902E-11 

3.772658 

1200 


1.27499E-U 

3.772658 

Plastic 

Strain = f (Temperature & Pressure) 


e p = A + BT + CP 

where. 


A = 03.1091E-3 
B = 3.96503E-7 
C = 2.96807E-3 


















Table VI. Comparison of the Plastic Strain Prediction Methods. 


Load Cases 

Synthesized 

Computed 

Percent 

Error 

• Example A - Nonlii 

1 

2 

3 

4 

lear Response Comput 

561 

892 

1356 

1377 

ed as a Function c 

576 

908 

1358 

1360 

>f Stress 

-2.6 

-1.8 

-0.1 

+1.2 

• Example B - Nonlii 

1 

2 

3 

4 

lear Response Comput 

576 

914 

1354 

1366 

ed as a Function c 

576 

908 

1358 

1360 

>f Temperature 
0 

+0.6 

-0.3 

+0.4 

• Example C - Nonlii 
Strai| 

1 

2 

3 

4 

lear Response Comput 
ght Linear Interpola 

576 

907 

1367 

1366 

ed Directly From I 
tion 

576 

908 

1358 

1360 

)ata Base Using 
0 

-0.1 

+0.7 

0 

• Example C - Nonlii 
Pressi 

1 

2 

3 

4 

lear Response Comput 
are 

562 

938 

1354 

1343 ! 

ed as a Function c 

576 

908 

1358 

1360 

>f Temperature/ 

-2.4 

+3.3 

-0.3 

-1.2 
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The shingle combustor segment model was generated having a thermal dis- 
continuity, and this hot spot is our critical location. In our current work, 
we have chosen to impose elastic constraints on the region of local inelas- 
ticity. The inelastic region is surrounded on all sides by purely elastic 
material, and it is our intention to examine overall structural and material 
behavior as the size of this critical location is altered. 

In our analysis, incremental cyclic plasticiity is performed using the 
Besseling subvolume method and employing the Von Mises yield criteria and the 
Prandtle Ruess flow rule. A strain-hardening creep law equation of the 
following form is used to describe the creep response of the material. 

e = ka n t m + qa r t 
c 

where the constants m, n, k, q and r are material and temperature dependent 
creep coefficients. 

It is our desire to determine if the size of the nonlinear region has 
significant effect on our stress-strain decomposition and synthesis tech- 
niques. Reduction in the number of nonlinear elements in the model can pro- 
duce significant savings in computational costs associated with a nonlinear 
analysis. For such an investigation the combustor shingle model was run for 
two separate analysis conditions. One being that the inelastic region is 
defined by eight elements and the other being that the inelastic region is 
defined by a single element. These two test conditions were each run three 
separate times to separate the effects of creep and plasticity. 

The shingle segment is shown in Figure 9 and modeled as illustrated in 
Figure 10. The thermal condition of the combustor shingle at peak temperature 
is shown in Figure 11(a) and the thermal cycle at the center of the hot spot 
is presented in Figure 11(b). The hot spot encompasses 12 elements with 1 
element (No. 99) exhibiting the maximum nonlinear material behavior (or total 
strain response) . 

The baseline condition for the current model was contained in NASA 
CP-2271. The baseline case represents the condition when all elements in the 
model are capable of nonlinear material behavior. Noting that we treated the 
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inelastic region as being one of eight elements in size and that the defined 
critical element denotes behavior of the model, we can observe the effects of 
material nonlinearity. Tables VII through XII give the results at the criti- 
cal element for the associated test conditions. 

Case 1: Baseline - All elements may be inelastic. 

Case 2: Inelastic region is eight elements in size: Elements 75-78, 

97-100 (See Figure 12) . 

Case 3: Inelastic region is one element in size: Element 99 (See 

Figure 12). 

It can be observed that when compared to the baseline case, original model 

behavior is best approximated when the inelastic region is eight elements in 

/ 

size. 

In support of the functional decomposition techniques, we investigated 
the influence of mechanical loads on our baseline models. The baseline models 
are those presented in NASA CR-165268 and NASA CP-2271. Tables XIII and XIV 
give pertinent data pertaining to the superposition of the mechanical loads. 
Simple pressure loads were sued for such purposes. 

In reference to Table XIII, the turbine blade tip model: 

D = Original model with imposed displacement boundary conditions. 

P = Original model with simple pressure load (8000 psi) . 

DP = Original model with imposed displacement boundary conditions 

and simple pressure load. 

In reference to Table XIV, the shingle combustor model: 

E = Elastic only 
ED = Elastic/Plastic 
EPC = Elastic/Plastic/Creep 

EPCP = Elastic/plastic/Creep With 1500 psi Pressure Load 

In support of the functional decomposition and synthesis techniques, we 
have investigated the influence of creep time on our baseline models. Six 
cases were examined and are presented here. The model used was the shingle 
combustor model and creep times were increased from 0.05 hour to 1.0 hour. 
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R-Total Strain 



i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

16 

EC 

-122 

-780 

-978 

-1184 

-1963 

-1113 

-460 

-975 

-656 

-633 

-380 

-809 

-558 

-544 

-321 

-314 

EP 

-123 

-835 

-1092 

-1420 

-2601 

-1888 

-1003 

-1474 

-1155 

-1131 

-879 

-1318 

-1070 

-1056 

-833 

-826 

EPC 

-123 

-835 

-1092 

-1420 

-2736 

-2037 

-1126 

-1582 

-1262 

-1239 

-986 

-1420 

-1171 

-1157 

-934 

-926 

R-Plastic 

Strain 















EC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

0 

-113 

-208 

-448 

-1516 

-1635 

-1053 

-1129 

-1129 

-1129 

-1129 

-1170 

-1170 

-1170 

-1170 

-1170 

EPC 

0 

-113 

-208 

-448 

-1516 

-1650 

-1099 

-1118 

-1118 

-1118 

-1118 

-1136 

-1136 

-1136 

-1136 

-1136 

R-Creep Strain 















EC 

0 

0 

0 

0 

-736 

-736 

-736 

-736 

-736 

-736 

-736 

-736 

- 736 

-736 

-736 

- 736 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

. 0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

0 

0 

0 

-301 

-301 

-301 

-301 

-301 

-301 

-301 

-301 

-301 

-301 

-301 

-301 




R-Total Strain 



i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

16 

EC 

-122 

-780 

-978 

-1184 

-1996 

-1122 

-463 

-982 

-663 

-640 

-388 

-819 

-568 

-554 

-331 

-327 

EP 

-122 

-841 

-1081 

-1375 

-2152 

-1325 

-469 

-990 

-672 

-649 

-396 

-834 

-583 

-569 

-346 

-342 

EPC 

-122 

-841 

-1081 

-1375 

-2272 

-1326 

-551 

-1069 

-750 

-727 

-475 

-906 

-656 

-642 

-418 

-414 

R-Plastic 

Strain 















EC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

0 

-120 

-209 

-402 

-1107 

,-1190 

-695 

-708 

-708 

-708 

-708 

-739 

-739 

-739 

-739 

-739 

EPC 

0 

-120 

-209 

-402 

-1107 

-879 

-618 

-618 

-618 

-618 

-618 

-625 

-625 

-625 

-625 

-625 

R-Creep Strain 















EC 

0 

0 

0 

0 

-759 

-759 

-759 

-759 

-759 

-759 

-759 

-759 

-759 

-759 

- 759 

-759 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

0 

0 

0 

-278 

-278 

-278 

-278 

-278 

-278 

-278 

-278 

-278 

-278 

-278 

-278 






R-Total Strain 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

EC 

-122 

-780 

-978 

-1184 

-2231 

-1369 

-712 

-1227 

-908 

-885 

-633 

-1061 

-811 

-797 

-574 

-567 

EP 

-122 

-866 

-1143 

-1517 

-2825 

-2121 

-1177 

-1727 

-1408 

-1385 

-1133 

-1578 

-1327 

-1313 

-1090 

-1083 

EPC 

-122 

-866 

-1143 

-1517 

-3088 

-2492 

-1473 

-1997 

-1680 

-1657 

-1405 

-1841 

-1591 

-1597 

-1354 

-1346 

R-Plastic 

Strain 















EC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

- 

“128 

-234 

-480 

-1702 

-1849 

-1268 

-1352 

-1352 

-1346 

-1352 

-1402 

-1402 

-1402 

-1402 

-1402 

EPC 

0 

-128 

-234 

-480 

-1702 

-1984 

-1399 

-1432 

-1432 

-1432 

-1432 

-1454 

-1454 

-1454 

-1454 

-1454 

R-Creep Strain 















EC 

0 

0 

0 

0 

-852 

-852 

-852 

-852 

-852 

-852 

-852 

-852 

-852 

-852 

-852 

-852 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

0 

0 

0 

-349 

-349 

-349 

-349 

-349 

-349 

-349 

-349 

-349 

-349 

-349 

-349 




Z -Total Strain 



i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

EC 

-835 

-1643 

-1885 

-2144 

-3301 

-2260 

-1465 

176 

-373 

- 459 

-938 

268 

410 

-224 

-618 

-613 

EP 

-834 

-1753 

-2115 

-2565 

-4084 

-3171 

-2341 

-551 

-1100 

-1187 

-1667 

-402 

-810 

-877 

-1272 

-1267 

EPC 

-834 

-1753 

-2115 

-2565 

-4389 

-3507 

-2527 

-824 

-1373 

-1460 

-1940 

-683 

-1092 

-1159 

-1554 

-1548 

Z-Plastic 

Strain 















EC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

0 

-325 

-549 

-1027 

-2635 

-2825 

-2366 

-2197 

-2197 

-2197 

-2197 

-2062 

-2062 

-2062 

-2062 

-2062 

EPC 

0 

-325 

-549 

-1027 

-2635 

-2783 

-2176 

-2072 

-2072 

-2072 

-2072 

-1960 

-1960 

-1960 

-1960 

-1960 

Z-Creep Strain 















EC 

0 

0 

0 

0 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

0 

0 

0 

-472 

-472 

-427 

-472 

-472 

-472 

-472 

-472 

-472 

- 472 

-472 

-472 




Z-Total Strain 



i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

EC 

-835 

-1643 

-1885 

-2144 

-3556 

-2506 

-1706 

-66 

-615 

-701 

-1181 

25 

-400 

-467 

-862 

-857 

EP 

-835 

-1817 

-2191 

-2652 

-4170 

-3249 

-2365 

-553 

-1101 

-1188 

-1667 

-374 

-799 

-865 

-1260 

-1255 

EPC 

-835 

-1817 

-2191 

-2652 

-4493 

-3566 

-2532 

-819 

-1359 

-1445 

1925 

- 642 

-1067 

-1133 

-1528 

-1522 

Z-Plastic 

Strain 


• 













EC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

0 

-357 

-601 

-1059 

-2753 

-2935 

-2446 

-2239 

-2239 

-2239 

-2239 

-2081 

-2081 

-2081 

-2081 

-2081 

EPC 

0 

-357 

-601 

-1059 

-2753 

-2823 

-2219 

-2082 

-2082 

-2082 

-2082 

-1947 

-1947 

-1947 

-1947 

-1947 

Z-Creep Strain 










• 





EC 

0 

0 

0 

0 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

-1562 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

0 

0 

0 

-472 

-472 

- 472 

-472 

-472 

-472 

-472 

-472 

- 472 

-472 

- 472 

-472 



























Z-Total Strain 



i 

2 

3 

4 

5 

6 

7 

8 

EC 

-835 

-1643 

-1885 

-2144 

-3301 

-2260 

-1465 

176 

EP 

-834 

-1753 

-2115 

-2565 

-4084 

-3171 

-2341 

-551 

EPC 

-834 

-1753 

-2115 

-2565 

-4389 

-3507 

-2527 

- -824 

EPCP 

-866 

-1787 

-2150 

-2603 

-4370 

-3512 

-2543 

-813 

Z-Plastic 

Strain 







EC 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

0 

-325 

- 549 

-1027 

-2635 

-2825 

-2366 

-2197 

EPC 

0 

-325 

'-549 

-1027 

-2635 

-2783 

-2176 

-2072 

EPCP 

0 

-328 

-557 

-1017 

-2584 

-2806 

-2190 

-2091 

Z-Creep Strain 







EC 

0 

0 

0 

0 

-1420 

-1420 

-1420 

-1420 

EP 

0 

0 

0 

0 

0 

0 

0 

0 

EPC 

0 

0 

0 

0 

- 472 

-472 

-472 

-472 

EPCP 

0 

0 

0 

0 

-482 

-482 

-482 

-482 


Z Components of Strain 


9 

10 

11 

12 

13 

14 

15 

16 

-373 

-459 

-938 

268 

410 

- 224 

-618 

-613 

-1100 

-1187 

-1667 

-402 

-810 

-877 

-1272 

-1267 

-1373 

-1460 

-1940 

-683 

-1092 

-1159 

-1554 

-1548 

-1391 

-1478 

-1958 

-700 

-1110 

-1177 

-1571 

-1563 

0 

0 

0 

0 

0 

0 

0 

0 

-2197 

-2197 

-2197 

-2062 

-2062 

-2062 

-2062 

-2062 

-2072 

-2072 

-2072 

-1960 

-1960 

-1960 

-1960 

-1960 

-2091 

-2091 

-2091 

-1980 

-1980 

-1980 

-1980 

-1980 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

-1420 

0 

0 

0 

0 

0 

0 

0 

0 


-4 72 -472 -472 -472 -472 -472 -472 -472 

-482 -482 -482 -482 -482 -482 -482 -482 






Table XIV. Case 1 Results - R Components of Strain (Concluded) 


R-Total Strain 



i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

EC 

-122 

-780 

-978 

-1184 

-1963 

-1113 

-460 

-975 

-656 

-633 

-380 

-809 

-558 

-544 

-321 

-314 

EP 

-123 

-835 

-1092 

-1420 

-2601 

-1888 

-1003 

-1474 

-1155 

-1131 

-879 

-1318 

-1070 

-1056 

-833 

-826 

EPC 

-123 

-835 

-1092 

-1420 

-2736 

-2037 

-1126 

-1582 

-1262 

-1239 

-986 

-1420 

-1171 

-1157 

-934 

-926 

EPCP 

12 

-680 

-921 

-1220 

-2377 

-1612 

-697 

-1289 

-851 

-828 

-576 

* 

-1019 

-769 

- 756 

-533 

-533 

R-Plastic 

Strain 















EC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EP 

0 

-113 

-208 

-448 

-1516 

-1635 

-1053 

-1129 

-1129 

-1129 

-1129 

-1170 

-1170 
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1. All elements in the model are capable of nonlinear material 
behavior. 

2. An isolated patch of elements is capable of nonlinear material 
behavior. 

3. A single element is capable of nonlinear material behavior. 

Shown in Tables XV and XVI are the results of the study. 

Case 1 - Elastic/Plastic/Creep Behavior 

Case 2 - Elastic/Creep Behavior Only 

We also examined the following two conditions: 

1. All elements in the model are capable of nonlinear material 
behavior. 

2. An isolated patch of eight elements are capable of nonlinear 
material behavior. The nonlinear region is surrounded by purely 
elastic material. 

For the above conditions the analysis was performed using creep times of 

0.05, 0.1, 0.5, and 1.0 hours. The results are summarized in Tables XVII 
through XX. 

The following outline summarizes our thoughts on the direction of the 
decompositon and sysnthesis methods. 


II. DECOMPOSITION TECHNIQUES 

A. Component Identification 

1. Will concern only the following components: 

• Combustor Liner 

• Turbine Blade 

• Turbine Vane 

2. No mesh generation or mesh refinement will be permitted once 
decomposition procedures have begun. 

B. Thermomechanical Loads 

1. Operating Conditions are identified by temperatures, pressures, and 
rotational speed. 


43 



R-Total Strain 


(a) 

1 

-123 

2 

-835 

3 

-1092 

4 

-1420 

5 

-3636 

6 

-2654 

7 

-1941 

8 

-2402 

9 

-2082 

10 

-2059 

11 

-1807 

12 

-2227 

13 

-1974 

14 

-1960 

15 

-1738 

16 

-1725 

(b) 

-122 

-866 

-1143 

-1517 

-4430 

-3244 

-2489 

-2970 

-2649 

-2626 

-2374 

-2783 

-2532 

-2518 

-2296 

-2881 

(c) 

-122 

-780 

-978 

-1183 

-2657 

-1789 

-1127 

-1637 

-1318 

-1295 

-1043 

-1469 

-1218 

-1204 

-981 

-971 

R-Creep Strain 















(a) 

0 

0 

0 

0 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

-1683 

(b) 

0 

0 

0 

0 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

-2204 

(c) 

0 

0 

0 

0 

-2286 

-2286 

-2286 

-2286 

-2286 

-2286 

-2286 

-2286 

-22 86 

-2286 

-2286 

-2286 

Z-Total Strain 















(a) 

-834 

-1753 

-2115 

-2565 

-5886 

-4527 

-3608 

-1951 

-2500 

-2587 

-3066 

-1817 

-2232 

-2299 

-2693 

-2681 

(b) 

-835 

-1817 

-2191 

-2652 

-5857 

-4476 

-3572 

-1877 

-2423 

-2509 

-2988 

-1724 

-2149 

-2215 

-2610 

-2598 

(c) 

-835 

-1642 

-1885 

-2144 

-3993 

-2958 

-2158 

-510 

-1059 

-1146 

-1625 

-416 

-841 

-968 

-1303 

-1297 

Z-Creep Strain 















(a) 

0 

0 

0 

0 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

-2666 

(b) 

0 

0 

0 

0 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

-2591 

(c) 

0 

0 

0 

0 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 

-3578 





R-Total Strain 


(a) 

1 

-122 

0 

~ g 

1 

3 

-978 

4 

-1184 

5 

-3215 

6 

-2416 

7 

-1774 

8 

-2276 

9 

-1957 

10 

-1934 

11 

-1682 

12 

-2102 

13 

-1851 

14 

-1837 

15 

-1615 

16 

-1597 

(b) 

-122 

-780 

-978 

-1184 

-4205 

-3397 

-2751 

-3247 

-2928 

-2905 

-2653 

-3071 

-2820 

-2806 

-2584 

-2565 

(c) 

-122 

-840 

-1081 

-1375 

-2679 

-16i9 

-910 

-1422 

-1104 

-1081 

-828 

-1249 

-998 

-983 

-761 

-754 

R-Creep Strain 















(a) 

0 

0 

0 

0 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

-2739 

(b) 

0 

0 

0 

0 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

-3632 

(c) 

0 

0 

0 

0 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 

-1233 


Z-Total Strain 


(a) 

-835 

-1643 

-1885 

-2144 

-5341 

-4370 

-3599 

-1947 

-2494 

-2580 

-3059 

-1847 

-2271 

-2338 

-2732 

-2717 

(b) 

-835 

-1643 

-1885 

-2144 

-5732 

-4756 

-3978 

-2323 

-2871 

-2957 

-3436 

-2223 

-2648 

-2714 

-3109 

-3095 

(c) 

-835 

-1762 

-2076 

-2467 

-4018 

-2720 

-1858 

-202 

-749 

-835 

-135 

-73 

-496 

-563 

-957 

-952 


Z-Creep Strain 


(a) 

0 

0 

0 

0 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

-4706 

(b) 

0 

0 

0 

0 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

-5159 

(c) 

0 

0 

0 

0 

-1643 

-1643 

-1643 

-1643 

-1643 

-1643 

-1643 

-1643 • 

-1643 

-1643 

-1643 

-1643 















































2. These operating conditions define the individual misison points for 
each component and are intended to encompass the expected range of 
engine operation. 

C. Component Analysis 

1. Total structural response of the generic component will be obtained 
through detailed finite element anlaysis procedures. 

2. A detailed analysis will be performed for each set of operating 
conditions (mission points) and critical locations can be identified 
for the component under study. This procedure can be user-defined or 
automated. 

3. Critical Locations - Once the critical locations are identified, 
procedures will be implemented to store information on (a) nodal 
displacements, and (b) elemental stresses and strains for the nodes 
and elements defining the region. This information will form the 
data base on which the synthesis techniques will be based. 

D. Structural Response and Load Application 

1. The structural response at the critical location can fall into these 
categories: 

• Uniaxial 

• Biaxial 

• Triaxial 

2. The generic mode of load application can be identified and a 
hierarchy defining the overall influence and effect on overall 
structural response can be formed. 

E . Material Response 

1. Material response can be modeled in a number of ways. The following 
limitations can be imposed on the material behavior: 

• Elastic only 

• Elastic/plastic 

• Elastic/creep 

• Elastic/plastic/creep 

• Flow rule and yield criteria 

• Creep formulation 

F. Decomposition Procedures and Models 

1. Serving as baseline data sets are models previously validated under 
other NASA programs. 
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• Turbine blade tip durability 

• Inelastic 

2. These models will be used to create a preliminary data base from 
which we will evaluate our stress/strain decomposiiton techniques. 

3. These models represent actual engine conditions and are intended to 
economically and efficiently evaluate the stress/strain 
decomposition and synthesis techniques. 


II. SYNTHESIS TECHNIQUES 


A. Component-Specific Synthesis Techniques 

1. The synthesis techniques will be based on a defined set of mission- 
segraent/component station characteristics. 

• The comnponent station characteristics are defined in the 
thermomechanical load (TML) matrix. 

• The structural data base will supply the necessary mission 
segment response upon which the synthesis alrogithms will 
operate. 

B. Synthesis Algorithms 

1. Algorithms used to compute fixed point mission-time profiles of 
"local" (critical point location) stresses, strains, and 
displacements. These algorithms will be limited to local level 
response. These algorithms may be based on empirical relations 
between loads and deformation states. 


III. COMPONENT DATA BASE STRUCTURE 
A. Component Dependent 

1. Three distinct data bases will arise from the component specific 
modeling effort: 

• Turbine blade data base (TBDB) 

• Turbine vane data base (TVDB) 

• Combustor data base (CBDB) 

2. The component data bases will be identical in structure but 
dissimilar in content. 

3. A simple set of access routines will be used to construct and access 
the three data bases. 
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Component Specific 
Modeling Effort 


T 1 

¥ 1 

¥ 1 

Turbine Blade 

Turbine Vane 

Combustor 

Data Base 

Data Base 

Data Base 


B. Component Data Base Attributes 

1. Random access 

2. Self-contained sizing 

3. Simple structure 

4. Three consitutent parts 

• Control includes 

- Overall problem sizing 

- Critical region and mission point specifications 
Assorted pointers and counters 

- Descriptions and titles 

• Critical Region Specification Includes 

Location and number of critical points are defined 

- Assorted pointers and counters for the mission points 
are defined. 

• Mission Point Specification Includes 

- Local control information 

- Node and element information and data 
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The following sketch shows how three constituent parts interact. 



General Data Base Detail: Mission Point 

Structure. Specification 


2.6 TASK II - STRUCTURAL ANALYSIS METHODS EVALUATION 

This task was completed during 1983. The selections outlined in 2.1 are 
the result of this effort. 
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2.7 TASK VII - THERMODYNAMIC LOADS MODEL 

The general form of the input, output, and calculations for blades, vanes, 
and combustor thermodynamic loads was defined during 1983, and work was initi- 
ated on the software coding. The final software specification is shown in the 
following paragraphs. 

Thermodynamic Loads Features for TDE Model: 

I. BLADES AND VANES 


A. Input will be : 

1. A series of midspan station temperatures at a specified SS ref 
erence case condition, and for each station a percent cord 
envelope dimension 

2. A set of output percent cord envelope dimensions 

3. A table of percent radial span versus cooling effectiveness 
factors. RF, where: 

ri 

“'jr 2 ■ 

'cms 

H cs = cooling effectiveness at specified span dim. 

^csra = co °li n 8 effectiveness at midspan 

The point density in this table will be such that spanwise 
linear interpolation will suffice. 

B. Calculations 

1. For each temperature a nominal cooling effectiveness will be 
calculated from: 


_ T4l-Tm 
c T41-T3 
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2. At each new condition, all station r| values will be modified 
by a factor, CF, as follows: 

a - n c ) 

^ ” ( i _ q )REF - ^ / (T3 ppp ,v T4l) ) a , where a is input 

3. Output n c values will be linearly interpolated based on percent 
cord envelope at adjacent input stations 

4. For each output radial distance (specified in the input) all 
station temperatures will be calculated as follows: 

Tij = T3 + (1 - n )*CF(T4l-T3)*RF 
C R 

Where i = station index 

j = radial station index 


II. COMBUSTOR 

A. Inputs will be : 

1 . Metal bulk made temperatures at axial stations on inner and 
outer liners, and for each station, x and y dimensions of node 
centers at a specified reference case condition. Both hot 
streak and average metal temperature values will be input. 

2. It is assumed that output node dimensions will match the input 
nodes. If cross-meshing is to be required, it will be done 
before the input is defined. 

3. A set of linear equation constants will be input for a sparse 
set of axial locations, identified by "AT station" numbers. 
Output AT locations will be same as input locations. 

4. A set of AP scaling constants will be inpout at a sparse set of 
locations, identified by AP station numbers. Output AP stations 
will be same as input. 

B. Calculations 


1. A nominal cooling effectiveness will be calculated for each 
input combustor average and hot streak temperature. 

2. At each flight condition, each metal temperature will be recal- 
culated from: 
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Tm. . = T3. + (1 - n )(T41.-T3.) 

ij J ^ J J 

where i = station index and j = flight phase index 
3. At each flight condition, each AT will be recalculated from: 

AT.. = [T3.+ (1 - n )(T4l-T3)][b.+m.*P3.] 
ij J 'c. i i 


4. At each flight condition, each AP will be recalculated from: 


AP. . = P3.*K.* 
J i 



The accuracy of the thermodynamic engine model has been evaluated, rela- 
tive to the engine steady-state performance computer model (cycle deck). In a 
model based on interpolation methods, the maximum error must occur in regions 
where the "distance” from known data points is greatest. Figure 13 shows the 
altitude versus Mach No. map of stored data points used in the thermodynamic 
engine model (TEM) . The X symbols on this map indicate the worst-case points 
selected for the error evaluation. 

The interpolation logic varies from quadrant to quadrant within a 
rectangle bounded by stored data points. The nearest point is always used as 
the base from which the interpolation process is started, for example. Figure 
14 shows a typical set of four stored data points. The shaded area is a 
quadrant. The evaluation test point is at the center. The hypothetical error 
magnitude curves drawn from two corners to the center illustrate two facts: 

• In any quadrant, the error surface is approximately parabolic in 
shape and maximum at the center of the stored data. 

• The four quadrants have different error surfaces, and a discontinuity 
occurs where they meet. 

Since there can be four different maximum error values at each test point, 
the error analysis was performed four times at each point. The results were 
summarized by a computer program. The right-hand four columns of Table XXI 
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Table XXX. Error Summary 
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AU6 


6.8828 

6.8859 

9.8878 

9.8689 

6.1125 

04 1664 

041119 

0.1199 






MAX 


-1.8963 

-1.8814 

-1.4891 

-1.3859 

-6.5570 

-9.6666 

-9.0497 

-6.9463 






HVAL 
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ai 

81 

15 

3 

3 

3 
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show the error values that exceed the target value which is listed in the 
center of the heading. The left-hand five columns identify the test points. 

The middle four columns show the accuracy level available before the 
improvements developed in this program were incorporated. The average of the 
absolute errors, the max error, and the number of "Exceedances" are given at 
the bottom of each column. Note that speed, pressure, and horsepower errors 
are expressed as percent of the rated standard day, sea level value. It seemed 
more meaningful to express temperature errors in degrees. 

Table XXII shows a brief summary of the accuracy level achieved. Column 2 
shows the average of all test point errors. Columns 3, 4, and 5 show the value 
that is exceeded 2, 4, and 11 times (1%, 2%, and 5% of the 220 error values). 
Note that all data in this figure refers to the worst-case test points. Since 
the error surfaces are approximately parabolic in shape, the average error in 
each quadrant is approximately half of the maximum error, and the overall error 
is approximately half of the average error listed. 


Table XXII. CF6-50C2 Validation Case 
Error Analysis. 



Error Exceeded N Times 

Average 

N = 2 

N = 4 

N = 11 

P 2 

0.03% 

<1% 

<1% 

<1% 

P 3 

0.23% 

1.3% 

1.2% 

1.0% 

FNIN1 

0.49% 

1.8% 

1.8% 

1.8% 

XN25 

0.17% 

0.7% 

0.6% 

0.5% 

T 

2 

0.08% 

<5° 

<5° 

<5° 

T 

3 

2.1° 

<10° 

<10° 

<10° 

*41 

11° 

47° 

47° 

35° 


The maximum error of 47° F listed for T4l may seem large. This was one 
of the most difficult parameters to fit. However, note that 47° is only 2% of 
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the rated T41 value (expressed in 0 F) and the true average error is only 
approximately 0.23%. 

Work is continuing toward the completion of the thermodynamic loads 
features to be incorporated in the thermodynamic engine model. 
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2.8 TASK VIII - COMPONENT SPECIFIC MODEL DEVELOPMENT 


2.8.1 Geometric Modeling 

A recipe for a segment of a combustor liner has been completed. Figure 
15 shows the typical nugget of a combustor liner. The recipe describes the 
geometry of a nugget as a function of a prescribed set of physical parameters. 
These parameters include the thicknesses of the liner skin, the radii of the 
internal and external surfaces of the nugget,, the height and length of the 
nugget, and the overall length of the panel. Generic master regions are then 
constructed using these key points. The coordinates of these key points and 
the connectivity of the master regions are then used as input into a 20-noded 
mesh generator to obtain the final mesh. In addition, the key information is 
written into an ESMOSS recipe format. This recipe is then used in conjunction 
with the ESMOSS mesh-generating logic to obtain a discretized mesh. Figure 16 
shows a typical nugget, some of the physical input parameters, and the master 
region definition based on these parameters. Figure 17 shows representative 
2D and 3D models generated from the master region. 

We came to the decision that we could not use exactly the same procedure 
for the turbine blade model. This is because of the significant differences 
in the two geometries. A typical turbine blade geometry is shown in Figure 
18. The outside of this airfoil is a complex curve defined by aero require- 
ments. There are no reasonable number of physical parameters which could be 
used to describe this geometry. We evaluated many options and believe that 
the best course is to use a data file of coordinates which define this closed 
curve. The capability to work with a file of points is present in ESMOSS. 

The interior of the cross section (wall thicknesses, location and number of 
spars, etc.) will be specified as physical parameters. The interior cross 
section can be defined using a data file or by specifying a constant wall 
thickness as a physical parameter. Figure 19 shows a plot of the exterior and 
interior points for a typical airfoil. 

Once the basic cross section is defined, the location and size of the 
interior spars can be controlled by the user through selection of paired 
points on the interior surface. The generator will then automatically insert 
the spars at the desired locations. The user can vary the number of spars 


60 




Typical Nugget 


Physical Input Parameters 



Master Region Definition 


Figure 16. Combustor Recipe Process. 
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2D Model 



3D Model 


Figure 17. Combustor Nugget Finite Element Models. 
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and the geometry of each spar. Figure 20 shows the pointer selected for a 
spar configuration. Figures 21 and 22 show coarse meshes for a blade with 
and without cavities. This software is complete and integrated with ESMOSS. 

2.8.2 Remeshing and Mesh Refinement 

This area and the next, Self-Adaptive Solution Strategies, touch on each 
other synergistically. What is sought in this program is the best combination 
of both. This involves two major areas of investigation: (1) the method to 

be used to refine, upgrade, and rearrange the mesh, and (2) the criteria to be 
used to activate this process. 

There are a number of ways to refine a mesh to get a better answer: 

(1) one way is to progressively subdivide a coarse mesh, always retaining all 
previous meshes within the finer mesh; (2) a second family of techniques 
totally realigns the mesh based on some criteria such as strain energy density; 
(3) a third method is to leave the mesh unchanged but upgrade the order of the 
elements . 

The first method, progressive subdivision, has certain theoretical and 
computational advantages. If the finite element interpolating functions used 
meet the requirements for completeness and continuity, convergence is mathe- 
matically guaranteed when we refine the mesh by progressive subdivision. The 
computational process of remeshing by progressive subdivision is straightfor- 
ward; however, it guarantees a larger problem to solve. 

For a solution of the finite element system of equations: 

[K] {6} = {F} 

suppose there is a numerical solution for the displacement, {6*}. Then the 
equilibrium or residual force vector is generated: 

{R} = {F} - [K] {w*} 

A perfect solution would result in this vector containing all zeros. 

Given the finite numerical accuracy of the computer, this is impossible. 
Therefore, a measure of the numerical "goodness" of the solution is to be 
found in how much this vector deviates from zero. Decisions on whether to 
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Figure 21. Coarse Mesh Selection for a Blade With Cavities. 
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Figure 22. Coarse Mesh Selection for a Blade Without Cavities. 

Ci 

to 



re-solve or redefine the problem can be based on the total and local devia- 
tions from zero. If a few local degrees of freedom are out of equilibrium, 
this might suggest a local remeshing. If the total equilibrium is deficient, 
this will require reraeshing and/or re-solving with greater numerical accuracy. 


The decision tree for this is as follows 


1 . 


If ZR. < C D 
1 R 

and all R. < CL. T 
1 RiL 


the solution is good 


2. If 2R. > CL 

i R 

and [Number of nodes with CR. T <R. < C p 1 > C 
then re-solve 1 1 iu s 

3. If IR. > CL 

l R 

and [Number of nodes with Ri > C p ] < C 
then remesh and re-solve K iu 15 

4. If R. < C D 

i R 

but some Ri > 

iu 

then remesh and re-solve 
where : 


R.^ = i fc k residual-free vector 

C n = Maximum allowable sum of R. 

R l 

= Lower bound for R. for possible remeshing 
iL 1 


C„ = Maximum allowable upper bound for R. 
iu 1 


Once an acceptable displacement solution has been reached, proceed to the 
element level. If, at the elastic level, stresses and strains are linearly 
connected, only one of these two needs to be evaluated. Strain will be 
checked. The total strain at each calculation point in an element is made up 
of an elastic strain and a thermal strain: 


e. = 

i 


1 


+ sT 


l 
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One aspect of this program is the establishment of acceptable strain 
gradients for different element types. Between adjacent strain calculation 
points in one element, and probably over the entire element, a strain gradient 
would not be chosen that could encompass an elastic-plastic-elastic or a plas- 
tic-elastic-plastic variation. Therefore, 

if 


e? - e® > | 2e . , , | , 

1 j - yield 1 ’ 

remesh this element. 

Additionally, there will be a change in sign, 
if 


Therefore, 



remesh this element. 

Once the nonlinear solution has been entered, the element level checks 
become more complex and more important. The total strain is now made up of 
the elastic strain, thermal strain, plastic strain, and creep strain: 


e. = e® + + eP + e 

11111 


Now stress and strain are no longer linearly connected; stress is a 
function of elastic strain only. Once again, between any two adjacent calcu- 
lation points within one element, an elastic strain gradient greater than the 
allowable material elastic gradient is not desirable. Thus, 

if 


£? - £ e . > | 2e . , , 

i j - yield 


remesh this element. The limit on the thermal strain would still be retained. 
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If 


— “ < 0, remesh this element. 

er 

J 

The next check is on the computed plastic and creep strain. No sign 
changes in either of these are allowed. In addition, a maximum gradient is 
set. 

If 


or 


or 


£? 

“ < 0 


— < 0 

e c 

J 


| £*?-£? 

1 .1 


I > c T 


or 


| £ ± - £j | > C c , 


remesh this element. 


Next, proceed to the interelement level check. These are of the same 
nature as the above, but now involve adjacent calculation points in adjacent 
elements. 

Our strategy for proceeding in this area consists of the following: 

• Select an analytical model with known solutions. 

• Use the "recipe procedure" to generate several meshes of different 
density. 

• Use the 20-noded finite element to evaluate remeshing criteria. 
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For our first model, we selected a semi-infinite plate with a circular hole 
near the edge. The stress concentration solutions for this problem for vari- 
ous loadings is known. Figure 23 shows the recipe parameters and master 
region definition for this problem. Figure 24 is a master region generated 
for a particular case. Figures 25 and 26 show a coarse discretized model 
based on this master region breakdown. Figures 27 and 28 show a refined dis- 
cretized model. With these models we will evaluate such correlative parame- 
ters as: 

• Strain energy density 

• Elemental stress/strain gradients 

• Interelement stress/strain comparison 

• Nodal stress comparison. 

In order to evaluate these criteria, three models were created for the 
problem of a hole with an edge effect in a semi-infinite plate. These three 
models represent three different mesh densities. The first load case run was 
a constant membrane stress condition applied to the right-hand edge of the 
model with symmetric conditions on the left-hand edge. This case has been 
run for all three mesh densities using 2D plane stress elements. Figures 29, 
30, and 31 are plots of the three models. 

The objective of this evaluation is to determine if the interelement 
stress distribution can be used to "monitor" the adequacy of the mesh density. 
An additional objective is to compare various methods of approximating the 
surface stresses based on the stress output at the centroid of the elements. 
These two efforts are synergistic since a good measure of the mesh adequacy is 
the correct surface stress and the surface stress is certainly a function of 
the element stress distribution. There are two approaches to approximating 
the surface stress from 2D elements. The first extrapolates the elemental 
stresses to the surface. The second method uses the nodal displacement at 
the surface to compute the surface strain. This surface stress is then com- 
puted based on this strain and the appropriate material properties. Figure 
32 shows the extrapolated stresses for each mesh density and the calculated 
stress based on strain. The vertical axis is stress and the horizontal axis 
is "wall thickness" which is the distance from the surface. 


73 



Master Region Model 


Parameters 


A = Height 

B = Relative Placement of Center 
of Hole Relative to A 
D = Diameter of Hole 

Input: D as a Percentage of the Height 

B as a Percentage of the Height 

Note: The Hole Cannot be Placed Beyond the 

Limits of the Beam 

(Ex) Input = 20%, 23.334% 

Set a 0.2A Hole 

0.2334 A Center 


Figure 23. Recipe Parameters and Master Region Model for Plate 
With Hole. 
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Figure 24. Generated Master Region. 
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Figure 26. Magnification of Hole Region for Coarse Model. 
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Figure 27. Refined Discretized Model. 
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Figure 28. Magnification of Hole Region of Refined Model. 
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Figure 29. Coarse Mesh. 
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Figure 30. Medium Mesh. 



Figure 31. Fine Mesh. 
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Figure 32. Effect of Mesh Density on Stress Prediction. 




The correct value obtained from Petersen is = 3.9. Note that the two 
coarse meshes underestimate the surface stress using both the extrapolation 
method and the surface strain method. The fine mesh model underestimates very 
slightly, but gets almost the same answer by both techniques. 

A fourth model was created and run for the tensile loading conditions. 

The basic geometrical model with the applied boundary conditions is shown in 
Figure 33. Also shown in Figure 33 are blow-ups of the mesh around the bottom 
of the hole for all fan mesh densities. 

The interelement stress distribution has been evaluated for these four 
models by approximating the stresses at the nodal points. These nodal point 
stresses were computed by weighted averaging of the elemental centroid 
stresses. The weighting coefficient is inversely proportional to the distance 
between the centroid and the node point. The maximum difference of the stress 
at any node was computed as the maximum difference of any of the centroid 
stresses of elements adjacent to the node. The percent difference is a meas- 
ure of the mismatch of nodal stress normalized by the local average stress. 
Figure 34 is a tabulation that represents the maximum difference and percent 
difference of nodal stress for seven points that are common in all four 
models. Notice that for all the points,, the maximum stress difference and 
the percent difference decreases as the mesh density gets finer. This 
indicates that the percent difference could be used as a mesh refinement 
criteria. Further work is needed to determine if this criteria holds for 
other loadings and different models and to select the acceptable levels of the 
percent difference at which no further refinement is required. 

Different loading conditions are being generated for these models. In 
addition, an 8-noded and 20-noded 3D model are also being created. 

2.8.3 Self-Adaptive Solution Strategies 

We have successfully developed and incorporated a dynamic time increment- 
ing scheme into a version of the two-dimensional nonlinear structural analysis 
program. The dynamic time incrementing technique is intended to provide for 
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convergence and promote stability in the nonlinear iteration schemes currently 
being used. 

In the analysis code, the constitutive model prepared by Bodner is being 
utilized to model time dependent inelastic flow. Bodner' s model requires that 
constitutive equations be numerically integrated over finite time intervals. 
This procedure is iterative as is the method of solution we use to solve the 
finite element equations. The dynamic time incrementing scheme developed is 
specific to these methods. 

In the implementation of the dynamic time incrementation scheme, the 
determination of an appropriate time increment can be computed three separate 
ways when the following simplified assumptions are made. 

1. Inelastic strain varies linearly with time - This allows the time 
increment to be computed by limiting the maximum change in the 
inelastic strain that can occur during the time interval. 


At i + 1 


= At. 


(Ae 1 ) n vi 
e allowable 


1 (Ae 1 ) 

e max 


where, 


At. _ = next time subincrement 
l+l 


At. = current time subincrement 
i 


(Ae ) - maximum inelastic strain increment occurring in the 

e max x ° 

current time step 


(Ae ) .. , , - maximum allowable inelastic strain increment, 

e aiiowabie T r , ^ 

User input or set to default value. 


2. Stress varies linearly with time - This allows the time increment to 
be computed by limiting the maximum change in stress that can occur 
during the time interval. 


(A6 ) n __ 

At.*, = At. , A .V U °“ able 

i+l i (A6 ) 


e max 
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where 


(AS ) 
e max 


= maximum change in stress occurring in the current 
time step 


(A ^ e ) a ;Q owa bi e = maximum allowable stress change. User input or set to 
default value. 


3. Control integration error - This allows for the control of the local 
integration error when computing the inelastic strain during the 
time increment. The local integration error is not allowed to 
exceed a predefined tolerance. 

it . f 12*™ + (^UogaMe l V3 

i+1 L T 1 J 

1 

where, 

TOL = user input 
• * *1 

£ = second derivative of the inelastic strain rate. 

The results in Table XXIII help illustrate the utility and stability of 
the dynamic time incrementing technique. In the generation of Table XXIII we 
illustrate the effect on creep response of Inco 100 when we control the maxi- 
mum allowable inelastic strain (As^) and the maximum allowable stress (Act) 
that can occur during any time increment. 

We have included in our nonlinear analysis code an alternate method for 
determining convergence of the numerical solution. A global Euclidian norm 
was defined as follows: 


** 12 
(eV 

hJ 

- 

J 

(e}-V 
-1 2 


N 1 2 
(07) 

J-l J 



< Tolerance 

N = No. of structural DOF 
i = Iteration No. 
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Table XXIII. Stability of Dynamic Time Incrementing. 









Material Constants 








E - 21.3E6 psi 

(a). 

Strain 

Control. 





D - 10 4 Sec-1 








N - .7 








Zo - 915 ksi 


0 

130 


130 


130 

Zi - 1015 ksi 


if t 

126 


628 


5751 

Z2 * 600 ksi 

if 

Iteration 

252 


1256 


11502 

A - 0.0019 Sec-1 


A a 

4000 


4000 


2000 

r * 2.66 


AE 1 

.001 


.0001 


.00001 



CPU/ SNUB 

14.37/4U38 

58.44/4U040 

504.6/4U260 


t 

Et 

E 1 

e t 

E 1 

E T 

E 1 


10 sec. 

6136 

33 

6136 

33 

6136 

33 


5 min. 

12059 

5956 

12060 

5956 

12060 

5957 


10 

16034 

9931 

16039 

9936 

16039 

9936 


15 

19958 

13855 

19964 

13861 

19964 

13861 


20 

23880 

17777 

23886 

17782 

23886 

17783 


25 

27802 

21699 

27808 

21704 

27808 

21705 


30 

31724 

25621 

31730 

25626 

31730 

25627 


35 

35646 

29543 

35651 

29548 

35652 

29548 


40 

39568 

33465 

39573 % 

33470 

39574 

33470 i 


45 

43490 

37386 

43495 

37392 

43496 

37392 


50 

47412 

41308 

47417 

41314 

47417 

41314 


55 

51333 

45230 

51339 

45235 

51339 

45236 


60 

55255 

49152 

55261 

49157 

55261 

49158 


65 

59177 

53074 

59182 

53079 

59183 

53080 


70 

63099 

56966 

63104 

57001 

63105 

57001 

(b) 

. Stress 

Control. 







a 

130 


130 


130 



ift 

117 


126 


152 


if Iteration 

234 


252 


304 



A a 

6000 


4000 


2000 



AE 1 

.001 


.001 


.001 



CPU/SNUB 

13.76/4U464 

14.37/4U938 

16.96/4U406 


t 

Ex 

E 1 

e t 

E 1 

e t 

£l 


10 sec. 

6138 

35 

6136 

33 

6135 

32 


5 min. 

12059 

5956 

12059 

5956 

12059 

5956 


10 

16034 

9931 

16034 

9931 

16034 

9931 


15 

19959 

13855 

19958 

13855 

19958 

13855 


20 

23881 

17777 

23880 

17777 

23880 

17777 


25 

27803 

21699 

27802 

21699 

27802 

21699 


30 

31724 

25621 

31724 

25621 

31724 

25621 


35 

35646 

29543 

35646 

29543 

35646 

29543 


40 

39568 

33465 

39568 

33465 

39568 

33464 


45 

43490 

37387 

43490 

37386 

43490 

37386 


50 

47412 

41309 

47412 

41308 

47411 

41308 


-55 

51334 

45230 

51333 

45230 

51333 

45230 


60 

55256 

49152 

55255 

49152 

55255 

49152 


65 

59177 

53074 

49177 

53074 

59177 

53074 


70 

63099 

56966 

63099 

56966 

63099 

56996 
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Convergence is assumed when the difference in the norm of the unknowns 
(displacements in this instance) is less than or equal to some predefined 
tolerance. This concept represents a global test of convergence rather than 
a local or point- by-point method as we currently use. 

The model presented for illustration of such a technique is presented in 
Figure 35. Linear springs are used to connect nodes along the crack face and 
theoretically should carry no load when the crack has opened. The displace- 
ments of several nodes on this face are followed and presented in the follow- 
ing tables. The true elastic solution for such a model is located in the first 
column of each table and represents the solution attained when an infinitesimal 
tolerance is prescribed for convergence. For this model all stresses remain 
in the elastic region .and no time-dependent material behavior is assumed. 

For the model shown, a study was conducted to determine the quality of 
solution obtained when either a global or a local convergence criteria was 
used. Tables XXIV and XXVI are the results of utilizing local criteria and 
Table XXV when the global criteria is used. Attention should be directed to 
the number of iterations for solution as well as the quality, of the solution. 

We are continuing the effort to incorporate such concepts into our finite 
element code and evaluate their applicability to self-adaptive solution strate- 
gies currently under development. Our current work has focused upon vector 
norms due to their direct applicability and the ease with which they can be 
computed. The definition of the vector norm is essential to the development 
of techniques in error analysis for the vector norm provides the means by 
which quantitative and comparative statements of size and distance can be made 
between associated vectors. The theory and development of vector and matrix 
norms can be found in most texts on numerical methods and only the necessary 
concepts are presented here. 

The norm of an "N" dimensional vector is written as ||x|| , and represents 
a single number. The norm is a function of all the elements in X and should 
not be confused with the definition of the vector (Euclidian) length |X| . 

The norm of X satisfies the following conditions: 
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Table XXV. Global Euclidian Criteria 



Table XXVI. Percentage Change Criteria 


Node It 

Elastic 

10% 

5% 

1% 

. 50 

0.9224 

0.9193 

0.9193 

0.9174 

61 

1.575 

1.5685 

1.5685 

1.5661 

71 

2.1557 

2.1429 

2.1429 

2.142 

80 

2.7001 

2.6778 

2.6778 

2.6814 

90 

3.219 

3.1868 

3.1868 

3.1951 

100 

3.719 

3.6768 

3.6768 

3.6901 

No. of Iterations 
to Convergence 

(8) 

(8) 

(9) 

All displacements (X1C 

r 3 > 




9 






























where N is the null vector 


2. 1 1 otX 1 1 = |ct| ||x|| for all scales a and 

vectors X 

||X+y|| ||x|| + [[ Y || for all vectors X and Y 

There can be defined infinitely many norms but the three most common 
norms are the J2 : norms II • If p for P = 1, 2, and «. 

II * II 1 “ |Xil + |X 2 | + |X 3 | + |X N | 

II 'll 2 ■ [l*l| 2+ |x 2 | 2 + |X 3 | 2 + I%| 2 ] 1/2 

IML " max { Pi . |X 2 |, |X 3 | |X N |j 

The above norms are special cases of the vector norm: 



To utilize the norms defined above the concept of relative and absolute 
error can be defined (1) for any two vectors, the absolute error can be 
written as: 

e a = llxjl - llx-.JI 


And the relative error written as: 

e R Pill - IP-ill 
IP-ill 

The expression for relative error tends to be far more applicable to our needs 
and its structure lends itself conveniently to the formulation of convergence 



criteria in iterative methods. When compared to a predefined tolerance, the 
relative error reflects the percentage change that has occurred. 

To simplify the task of evaluating new convergence criteria, we have 
chosen to limit the present study to a class of problems that exhibit only 
geometric nonlinearity. We will use the nodal displacement as the control 
vector for which we compute the assorted vector norms, assemble the error 
expressions, and hopefully determine convergence. Noting that our solution 
scheme is an iterative technique, we are looking for the convergence of a 
series of vectors. For our application it is necessary to determine when a 
solution has been obtained as well as to determine the quality of the computed 
solution. 


Convergence of the numerical solution can be determined on either the 
local or global level. Currently within our code we use a local or point- by- 
point method and the present work represents an effort to see if a global 
scheme has applicability. Six separate convergence criteria were generated 
from the error expressions mentioned above and represent those most commonly 
used in the literature. 


Kill. - K-Ji 

1- Hii-i Hi 

PJ2 - \\x ± -i \\ 2 

3 - 


5. 


PiL - IIVilL 

Pi-ill. 


_ iix^-ikji, 

2 - Kill. 

P ± II 2 " Pi-lP 
4 - — rnm 

*1 2 


6 . 


KIL - K-ill. 

11*1 n„ 


In the above expressions, the subscripts are indicative of the iteration 

used. 

The model presented for the illustration of global convergence criteria 
is presented in Figure 35. This model is that presented in the March and 
April narratives. Linear springs are used to connect the nodes along the 
crack face and theoretically should carry no load when the crack has opened. 
We have limited ourselves to small strain/small displacements to reduce our 
test matrix. 
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The results of our work are presented in Tables XXVII, XXVIII, and XXIX. 
During each iteration the assorted vector norms are computed and the necessary 
error terms are formed. In all instances the solution has terminated when a 
local convergence criteria had been satisfied. The local criteria states that 
convergence can be assumed when during any iteration less than 0.5% change has 
occurred in the local displacement. As can be noted in the tables, the mea- 
sure of relative error is in excellent agreement with the measure of local 
convergence and thus the expression of relative error has potential applica- 
tion in our work as a measure of convergence when small displacement/small 
strain problems are encountered. 
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Table XXVII. Utilization of Vector Norms. 


Iterations 

Absolute Error 

Relative 

Error 

l|V ± II - II V-j-l II 

Pi || - ||v 1=1 || 

II Vi II - l|v 1= i 

| 

Vi-1 | 

l|v ± 

s 


f.Z6FZ!E W 

t.mm m 


f.318389E-f2 

P.:6t*4lE f? 

I.34129K ff 

5 

I.Z49334E-4Z 

f.v.mi ti 

f .247334E ff 

4 

I.193129E-I2 

». 7W1IE-H 

9.299931C ff 

7 

l.i32432E-f2 

f.5M722E-«: 

4.163294E ff 

* 

g 

f.U9fl9E-f2 

1.422954;-*: 

I.12739K ff 

V 

ft 

f.?2?ff7E-f3 

9.31M3K-K 

f.*52f2E-fl 

9 

If 

f.7Z5fl7E-f3 

I.239715H! 

f.774477E-fl 

11 

f.36373fM3 

I.192WM1 

l.6#6#4IC~f 1 

ft 

f.4414f3E~f3 

MW39M1 

I.472853E-II 

13 

f.3M368E-f3 

f.!r7Z4S-f! 

f.3689f6E~fl 

•• 

14 

f.268678E-#3 

I.S3;4flE*l2 

I.287822M1 

13 

f.2f9389E*f3 

f.:*-:37E-»2 

f.224S23E-fl 

14 

f.!633f9E-f3 

M«4?3E -ft 

f.l731S?E*fl 


Table XXVIII. Utilization of L(2) Vector Norms. 


Iterations 

Absolute Error 

Relative 

Error 

P ± II - llVi II 

II Vi II - Pi-1 || 

Pi II " Pi-i II 

ll*i-i II 

m ii 

t 

f.3W8fE-fl 

1.28*9*4; H 

9.382371E If 

4 

f.315477E-fl 

I.1725Z4E ff 

I.3I7744E ff 

3 

f.247281E~fl 

P.il5255E W 

f*24122«E ff 

4 

f.l92872E*fl 

f.9#8?$4E-fl 

9.188144E ff 

7 

f.l51f37E~fl 

I.584B2E-P1 

9.147334E ff 

1 

f.U8487E-fl 

§,43?363£*l: 

f.U3383E ff 

4 

•.424242E-IZ 

».3Z2*i9E-»l 

f.fft386E-fl 

If 

f.722Z71E-f2 

l.243744£-»l 

f.7f4544E*fl 

11 

f.364ff7E-«2 

f.l856KE-fi 

f.53fl82E-fl 

12 

f.43993fE*f2 

•.1*2233E-»I 

f.429163E~fl 

13 

«.34317fE-f2 

».if«i:E-f! 

f.334738E-fl 

14 

fl.2477«4E-«2 

f.84*3 s, .E-f2 

f.Z6U44E-fl 

13 

f.2f88f8E~f2 

•,633! 7 3E-I2 

f.2f3*9fE~fi 

16 

I.162887E-I2 

I.5KK1HZ 

f.!38894E-fl 
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Absolute Error 

Iterations || - Hv^ || 


3 

4 

5 
4 
7 

e 

9 

!# 

If 

12 

13 

14 

15 
1 * 


#.7115771-13 

I.354693M3 

#.432763E-#3 

#.337589E*f3 

•.2W337E-I3 

#.2#5415E-#3 

#.U#228E-«3 

I.124988E-I3 

#.9749i3E-#4 

l.76#492E-#4 

#.J932ME*#4 

9.442757E-I4 

#.36#952E-#4 

f.2815WE-#4 


of Vector Norms. 


Relative 

Error 

llv, II - 117,., II 

II Vi || - || Vi., || 

Vl-1 

llVi || 

I.2777I4C n 

i.lSRIIC* 

I.169647E M 

f.336#71E 91 

#.U3t7IE 94 

I.212199E «# 

f.793#ME-»l 

#.2f4534€ f# 

I.573175E-41 

9.159548E f# 

M228ME-91 

•.124455E ff 

#.31Wfc2E-#l 

#.97#77IMl 

I.239268E-91 

I.757264E*#! 

1. 182293E-91 

•.59W99E-91 

#.139**"E-#1 

f.4M758E-91 

#.1I7427E-#1 

#.3594#1E-I1 

#.829138E-42 

I.2M37IE-I1 

4.M1412E-I2 

9.218489E-I1 

M97154E-I2 

#.17#592E-I1 




APPENDIX A 


TASK II - DESIGN OF STRUCTURAL ANALYSIS 
SOFTWARE ARCHITECTURE 



u/;n/?s 


MODEL INO 
CRITERIA 


DECISION 

PYROMETERS 


CVCLE MAP 

MISSION DEFINITION 
ENGINE RATING 


PORT GEOMETRY 


HEAT XFER DATA 


CRITICAL POINT LOC 
STRUCTURAL DATA 


MATERIAL LIFE DATA 


MODEL 
SPEC I F IC 
COMPONENTS 


CRITICAL POINT 

iTfT 


CO 

^1 


COSMOS - CONTEXT 


0-0 


REVISION- 4 


8-10-03 



1 



Model Specific Component 








fwaa'21/79 



n't 


GENERATE PART HISTORY 


Kl VISION- 1 


7-22-0'J 












A4 . 1 


SET ANAL VS I S CONDITIONS 


REVISION- l 


O-10-B3 
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Analysis Program 


* 







PREPROCESSOR ATTRIBUTES: 


A. SOFTWARE 

• Modular in Structure 

• , Machine Independent 

• Low Core Requirement 

• Extensive Documentation 

B. FEATURES/OPERATION 

• User Friendly 

• Extensive Diagnostics 

• Complete Data Summary 

• Interactive Graphics 

- Model Geometry 

- Material Physical Properties 
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Mala 

Program 


Initialize 

Program 

Vectors 


Problem 

Definition 


Read & Store 
Mode Data 


Read 4 Store 
Element Data 


r Ust\ 
Load y 
.Case / Y «* 


Load Case 
Control 
Information 


Read 4 Store 
Material 
Physical 
Properties 


Read & Store 
Applied Nodal 
Displacement 


Read 4 Store 
Applied Nodal 
Forces 


Read 4 Store 
Nodal 

Temperatures 


Read 4 Store 
Element 

Pressure Loads 


Element 

Addition/ 

Removal 

Check 


Preprocessing 

Options 



Contour 

Plots 






POSTPROCESSOR ATTRIBUTES: 


A. SOFTWARE 

• Modular in Structure 

• Machine Independent 

• Low Core Requirement 

• Extensive Documentation 

B. FEATURES /OPERATION 

• User Friendly 

• Will Function in Batch or Time-Share Environment 

• Will Have Extensive Graphics Capabilities 

- X-Y Graph Plots 

Contour and Deflected Shape Plots at Defined Planes 

• Will Have Built-in Data Manipulation Routines (User Controlled) 

- User Friendly Data Base 






















APPENDIX B 


COMPONENT TEMPERATURE AND PRESSURE DECOMPOSITION 
AND SYNTHESIS PLAN 


j I. Blades and Vanes 


% 


A. Input will be: 


1. A series of mid-span station temperatures at a specified SS 
reference case condition, and for each station a percent cord 
envelope dimension. 

2. A set of output percent cord envelope dimensions. 

3. A table of percent radial span (or radii) versus temperature 
factors, RF, where: 


RF = 


Ts-T3 

Tms-T2 


Ts = temperature at specified span dim. 
Tms = temperature at midspan. 


The point density in this table will be such that spanwise linear interpo- 
lation will suffice. 


B. Calculations 


* 


1. For each temperature a nominal cooling effectiveness will be 
calculated from: 

_ T4l-Tm 
n c T41-T3 


2. At each new condition, all station r| values will be modified 
by a factor, CF, as follows: c 


CF 


(!-n c ) 

(l-n c )REF 


= (T3*T41)“/ (13^41)^ 


where a is input 


3. Output Nc values will be linearly interpolated based on percent 
cord envelope at adjacent input stations. 

4. For each output radial distance (specified in the input) all 
station temperatures will be calculated as follows: 
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Tij - T3 + (l-n CR )*CF(T41-T3)*RF 


Where i = station index 

j = radial station index 


II. Combustor 


A. Inputs will be : 

1. Metal bulk made temperatures at axial stations on inner and 
outer liners, and for each station, x and y dimensions of 
node centers at a specified reference case condition. Both 
hot streak and average metal temperature values will be input. 

2. It is assumed that output node dimensions will match the input 
nodes. If cross-meshing is to be required, it will be done 
before the input is defined. 

3. A set of linear equation constants will be input for a sparse 

set of axial locations, identified by AT station: numbers. 

Output AT locations will be same as input locations. 

4. A set of AP scaling constants will be input at a sparse set of 
locations, identified by AP station numbers. Output AP sta- 
tions will be same as input. 

B. Calculations 

1. A nominal cooling effectiveness will be calculated for each input 

combustor average and hot streak temperature. 

2. At each flight condition, each metal temperature will be recalcu- 
lated from: 

Tm.. = T3. + (l-n ) (T41.-T3.) 
ij J c i v j y 

Where i ■ station index 

j * flight phase index 

3. At each flight condition, each AT will be recalculated from: 



4. At each flight condition, each AP will be recalculated from: 

W41 2 

AP. . = P3.*K.* ~ 7 * T3. 
ij j l P3 j j 
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